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To be able to achieve their physics goals, future neutrino-oscillation experiments will need to 
reconstruct the neutrino energy with very high accuracy. In this work, we analyze how the energy 
reconstruction may be affected by realistic detection capabilities, such as energy resolutions, effi¬ 
ciencies, and thresholds. This allows us to estimate how well the detector performance needs to be 
determined a priori in order to avoid a sizable bias in the measurement of the relevant oscillation 
parameters. We compare the kinematic and calorimetric methods of energy reconstruction in the 
context of two —>■ disappearance experiments operating in different energy regimes. For the 
calorimetric reconstruction method, we find that the detector performance has to be estimated with 
an 0(10%) accuracy to avoid a significant bias in the extracted oscillation parameters. On the 
other hand, in the case of kinematic energy reconstruction, we observe that the results exhibit less 
sensitivity to an overestimation of the detector capabilities. 

PACS numbers: 14.60.Pq, 14.60.Lm, 13.15.-|-g, 25.30.Pt 


I. INTRODUCTION 

Long-baseline experimental searches of neutrino oscil¬ 
lations largely rely on the capability of pinning down the 
energy dependence of the oscillation probability, which is 
a nontrivial function of the true neutrino energy, E^,. As 
a consequence, the procedure employed to reconstruct 
the unknown incoming-neutrino energy from the mea¬ 
sured kinematics of the interaction products is a central 
element of the oscillation analysis. 

Experiments using neutrino beams peaked at E^, ^ 
600-800 MeV, such as T2K [T] and MiniBooNE [2], de¬ 
termine the energy distribution of charged-current (CC) 
events from the kinematics of the outgoing charged 
lepton—i.e., its kinetic energy and emission angle— 
measured by large Cherenkov detectors filled with wa¬ 
ter or mineral oil. This technique is mostly applied to 
quasielastic (QE) events—identified by the absence of pi- 
ons in the final state—that provide the dominant contri¬ 
bution to the total cross section at these energies. How¬ 
ever, it necessarily involves hypotheses on the reaction 
mechanism. 

The kinematic method of energy reconstruction is 
based on the assumptions that the beam particle inter¬ 
acts with a single nucleon at rest, bound with constant 
energy, and that no other nucleons are knocked out from 
the nucleus. It has been long known, however, that pro¬ 
cesses involving two-nucleon currents, final-state interac¬ 
tions, and nucleon-nucleon correlations give rise to the 
appearance of more complex final states, featuring more 
than one nucleon excited to the continuum. Reconstruc- 
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tion of the neutrino energy of such events, dubbed QE- 
like in the mesonless case [3], in general requires more 
complex methods, involving realistic models of nuclear 
dynamics. 

At energies well above 1 GeV, the contribution of in¬ 
elastic processes—resonant pion production and deep- 
inelastic scattering (DIS)—becomes larger and eventu¬ 
ally dominant. As a consequence, determination of the 
neutrino energy in this kinematic region requires recon¬ 
struction of events with many hadrons in the final state. 

As alternative to Cherenkov detectors in this regime, 
calorimeters measuring the visible energy associated with 
each event—i.e., the energy deposited by the final-state 
particles—have been proposed as effective devices, al¬ 
lowing for an accurate neutrino-energy reconstruction. 
Calorimeters are presently being used in the MINOS [3] 
and NOi^A [5] experiments. In their detectors, the total 
energy deposited by all reaction products is measured 
without a prior reconstruction of each final-state parti¬ 
cle’s track, momentum, or energy. The energy response 
for potentially complex final states is calibrated using 
test-beam exposures. 

On the other hand, future appearance experiments like 
DUNE [0] are expected to employ detectors capable of 
fine-grained tracking of a large number of interaction 
products. In their case, the tracking capability is the key 
to being able to select electron-neutrino events and dis¬ 
tinguish them from backgrounds, even for non-QE events. 

In this article, we compare the kinematic and calori¬ 
metric reconstruction methods in the oscillation analysis 
of a disappearance experiment. We also aim to explore 
the capabilities and limitations of a calorimetric analysis 
based on individually identified particle tracks and to de¬ 
termine what level of understanding of detector response 
and underlying events is required to meet certain physics 
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goals. 

For a muon-neutrino event, it is clear that the long 
muon track in itself is a clear signature, and no tracking 
beyond the leading muon is required for backgrounds re¬ 
moval. However, we study the disappearance channel 
because of its simplicity in terms of oscillation physics, 
and as a sandbox to develop suitable analysis tools. 

The calorimetric technique, as defined above, obviously 
rests on the ability of fully reconstructing the final state, 
which largely depends on the detector design and perfor¬ 
mance. Nuclear effects also play a role, as they may lead 
to a sizeable amount of missing energy, hindering the 
reconstruction of E^. For example, if a pion produced 
at the elementary interaction vertex is absorbed in the 
spectator system, in general, its energy is not deposited 
in the calorimeter. 

We consider an idealized setting, in which the near and 
far detectors are functionally identical, and a simple ex¬ 
trapolation between them can be performed [7] . To min¬ 
imize the uncertainty arising from nuclear interactions, 
the target nucleus selected for both detectors is carbon, 
the cross section of which has been extensively measured 
in a number of different channels 

The analyzed events are generated using the simulation 
code GENIE [29] supplemented with the vT package of ad¬ 
ditional modules [30], allowing us to describe the carbon 
ground state using the realistic spectral function m- 

To cover different experimental configurations, we con¬ 
sider two options: a low-energy (LE) setup, with a 
narrow-band off-axis beam, and a high-energy (HE) 
setup, featuring a broadband on-axis beam. In the LE 
option, the neutrino flux is peaked around 600 MeV, and 
the distance to the far detector is set to L = 295 km, 
while in the HE option, the flux is peaked at ^1-2 GeV 
and L = 1000 km. 

Event reconstruction is carried out assuming: 

(i) a perfect scenario, in which all produced particles 
are detected and their true energies are measured, 

(ii) a realistic one, in which detection efficiencies and 
thresholds are taken into account and the finite de¬ 
tector resolution leads to a smearing of the mea¬ 
sured energies. 

In both cases, neutrons are assumed to escape detection 
altogether. 

A key element of our analysis is the migration matrices, 
the elements of which, correspond to the probability 

that an event with a true neutrino energy in the jth bin 
ends up being reconstructed in an energy bin i. 

The event distribution representing the data is in each 
case generated using the realistic scenario. Fitted rates, 
on the other hand, are obtained using a linear combina¬ 
tion of migration matrices corresponding to the perfect 
and realistic scenarios. This procedure is very effective 
from the computational point of view, and provides a 
measure of the impact of detector performance on the 
oscillation-parameters fit. 


This article is structured as follows. In Sec. El we 
derive different methods of neutrino-energy reconstruc¬ 
tion, placing a special emphasis on the approximations 
involved. The neutrino cross section model employed for 
event generation and the treatment of detector effects 
are outlined in Secs. m and |IV[ respectively. In Sec. |V) 
we discuss the calculated migration matrices. Our oscil¬ 
lation analysis and its results are presented in Sec. lYll 
Finally, in Sec. IVII] we summarize our findings. 

II. ENERGY RECONSTRUCTION METHODS 

Consider CC neutrino scattering off a nuclear target, 
resulting in the knockout of n nucleons and production 
of m mesons. The energy and momentum conservation 
can be cast in the form 

Eu + Ma = Ei + Ea-u + E^p'+E^h', (1) 

^ 3 

k, = k^-p+Ep*+E^’ (2) 

i 3 

respectively, where E^ and {E^ and k^) are the neu¬ 
trino’s (charged lepton’s) energy and momentum, Epi 
and p( denote the energy and momentum of the ith 
knocked-out nucleon (1 < i < n), and E^' and h( stand 
for the energy and momentum of the jth produced meson 
(1 < j < m). The energy of the residual {A — n)-nucleon 
system, Ea-u, can be conveniently expressed as 

Ea-u = Ma — nM E + Ta-u, (3) 

in terms of the nucleon (target-nucleus) mass M (Ma), 
the recoil energy Ta-uj and the excitation energy E. In 
Eq. ([^, the recoil momentum of the system is denoted as 
— p, to allow the interpretation of p as the vector sum of 
the initial momenta of the knocked-out nucleons, assum¬ 
ing that it is not altered by final-state interactions with 
the residual system. 

Substitution of Eq. (§ into Eq. 0 leads to the neu¬ 
trino energy in the form 

E, = Ei + E + Ta-u + E(^P' - ^) + E 

i 3 

Note that while for mesons the total energies enter the 
sum, for nucleons only the kinetic energies contribute. 
This difference is a consequence of the fact that mesons 
are produced in the interaction process, whereas nucleons 
are only knocked out from the target nucleus. 

Assuming that multinucleon effects do not introduce 
strong energy dependence to the cross sections, the factor 
Cn = E + Ta-u can be treated as a constant at neutrino 
energies above several hundred MeV. Then, reconstruc¬ 
tion of the neutrino energy reduces to determining the 
energies of the particles in the final state |5^ . 

= Cn + Ei + E/('®p1 ~ ^ E/ 

i 3 
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This, so-called, calorimetric method can, in principle, 
be applied to any type of CC interaction. However, one 
needs to keep in mind that an accurate reconstruction of 
hadrons poses a formidable experimental challenge. In 
particular, neutrons typically escape detection, and any 
undetected meson results in energy underestimation by 
at least the value of the pion mass, ~135 MeV. 

When the invariant hadronic mass squared, defined as 

= (E^p; + E®-;)" - (Ep; + E^)’’ 


is known, Eqs. Q and ([^ can be solved for the neutrino 
energy, yielding the alternative expression [33] 


E„ = 


2EpEi - 2p ■ k, + W" - (El - p") - mj 
2{Ep — IpI cos 07 V — Ei + jk^l cos 9) 


( 6 ) 


where Ep = MA — EA-m 'me is the charged lepton’s mass, 
IpI COS 07V = P • ^v/Ey, and |k^| cos0 = k^ • \<i^/E^. 

On the other hand, in the case of a single-nucleon 
knockout associated with the production of m mesons, 
the requirement that the nucleon is on the mass shell, 
El, = + p'2, can be used to obtain the neutrino 

energy as 


E,, = 


M^ + jp-ke- E, - {Ep -Ee- E, E^f 
2[Ep - IpI cos 9n - Ee + |k^| cos 9 - Y,. Uj] 


(7) 

with T-Lj = Eh' - h' -k^j/E,,. 

In practice, an application of the above formulas re¬ 
quires (i) neglecting the unmeasured recoil momentum 
IpI and (ii) approximating the energy of the residual nu¬ 
clear system by a constant, which amounts to setting 
Ep = nM — €n in Eq. ([^ and Ep = M — e in Eq. Q. 
These simplifications lead to the expressions 


= 


2{nM — en)Ee + — {nM — — ' 

2{nM — Cn — Ee + |k^| cos0) 


( 8 ) 


and 


ET = 




(k, + ^ . h;)2 - (M - e - E, - E, Ee,, 


2{M-e-Ee + |k,| cos9 - ^ 


(9) 

Owing to difficulties with an accurate determination 
of the invariant hadronic mass, the use of Eq. (|^ is usu¬ 
ally restricted to the process of a single-nucleon knockout 
with no pions produced, in which is known and equals 
m 2 piiKllllMlISS]. On the other hand, Eq. has 
been applied in the energy reconstruction for single-pion 
events by the MiniBooNE Collaboration [12] , setting the 
single-nucleon separation energy e to zero. 

The kinematic energy reconstruction, by means of 
Eq. 0 or does not require the knocked-out nucleon’s 
momentum to be measured. However, as this method 
assumes specihc final states, its accuracy is spoiled by 
any undetected hadron. For example, when a produced 


pion is absorbed or undetected, the energy reconstructed 
from Eq. ^ under the QE hypothesis is typically lower 
than the true one by '^300-350 MeV; see Figs. 6 and 7 
of Ref. [Sni ■ While the process of multinucleon knockout 
affects the energy reconstruction in a similar way—i.e., it 
redistributes the strength of the reconstructed flux from 
the peak mainly to the low-energy tail ISZIISH] —this effect 
seems to be less relevant at higher beam energies |3S]. 

For completeness, we mention that neutrino energy can 
be also found exploiting momentum conservation only. 
Equation (§ multiplied by a factor k^jEy, 

Ej, = |k^| COS0 —IpI cos9n+''^ |Pi| cos0i-l-^ |h'| cos9j, 

i 3 

( 10 ) 

where |p'| cos9i = p'-kj,/E^ and |h'| cos0j = h' •kj./E^, 
shows that E^ can be determined from the projections of 
the final momenta on the beam direction. Neglecting the 
contribution of the recoil momentum, one obtains the 
expression 

E“°“ = |k^|cos0 + ^ Ip'I COS0* |h' |cos0j, (II) 

i 3 

a special case of which has been employed in an anal¬ 
ysis of single-nucleon knockout events by the NOMAD 
Collaboration [S|. 

Performing the kinematic energy reconstruction in this 
article, we employ Eq. Q assuming single-nucleon knock¬ 
out, regardless of the actual number of nucleons in the 
final state. We set to M^ for mesonless events and 
to M^, Ma = 1.232 GeV being the A resonance mass, 
when at least one meson is observed. The single-nucleon 
separation energy is fixed to 34 MeV. The same value of 
e is added in the calorimetric energy reconstruction (§ 
for every nucleon detected. 


III. EVENT GENERATION 

Our analysis is based on the description of nuclear 
structure and interaction dynamics in the genie Monte 
Carlo generator [21], version 2.8.0, supplemented with 
the lyT package of additional modules m- 

GENIE is a modern and versatile platform for neutrino 
event simulation. It has been developed putting special 
emphasis on scattering in the energy region of few GeV— 
important for ongoing and future oscillation studies— 
where various mechanism of interaction are relevant. A 
number of neutrino experiments employs this generator 
in data analysis HDj. 

Resonant pion production in genie, considered for 
IT < 1.7 GeV, is accounted for using the model of Rein 
and Sehgal m- Gompared to 18 in the original model, 
16 resonances of unambiguous existence are implemented 
using up-to-date parameters but neglecting the interfer¬ 
ence between them. The effect of the charged lepton’s 
mass is taken into account only in the calculations of the 
phase-space boundaries. 
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FIG. 1. (color online). Panel (a): Per-nucleon CC inclusive cross section divided by neutrino energy, obtained using genie 
2.8.0+ |29[I30) with the relativistic Fermi gas model (dashed line) and the spectral function approach (solid line) as a nuclear 

model in QE interaction. The results for carbon are compared to the experimental data for carbon extracted from NOMAD [S] 
and those for hydrocarbon (CH) reported from SciBooNE |14) and T2K [171119j . Panel (b): Breakup of the contributions to 
the inclusive cross section. The labels DIS, res, 2p2h, and QE refer to deep-inelastic scattering, resonant pion production, 
two-nucleon knockout, and quasielastic scattering, respectively. 


The contribution of nonresonant processes, classified 
in GENIE as DIS, is calculated following the method of 
Bodek and Yang [321113] • This effective approach extends 
the range of applicability of the parton model to low neu¬ 
trino energy by modifying the (leading-order) parton- 
distribution functions in the low-Q^ region, being 
the four-momentum transfer squared. Higher-order and 
target-mass corrections are accounted for by replacing 
Bjorken x with a new scaling variable (121 113] ■ While 
DIS in GENIE is the only mechanism of interaction at 
W > 1.7 GeV, it also produces one- and two-pion events 
in the resonance region. 

Hadronization in genie is performed using the 
Andreopoulos-Gallagher-Kehayias-Yang approach jH] , 
which combines Koba-Nielsen-Olesen (KNO) scaling [IS] 
at low values of the hadronic invariant mass with the 
pythia/jetset calculations [IS] at high W, ensuring a 
smooth transition between the two regimes. 

Two-nucleon knockout {2p2h) events are simulated in 
GENIE using the procedure of Dytman m, obtained 
modifying and extending the one of Ref. |3H] ? derived for 
electron scattering. The invariant mass of two-nucleon 
events is assumed to have a Gaussian distribution cen¬ 
tered a,t W = {M + Ma)/ 2. The charged lepton’s kine¬ 
matics is distributed according to the magnetic contri¬ 
bution to the cross section for scattering on a free nu¬ 
cleon |49j . identical for neutrinos and antineutrinos. The 
strength, adjusted to fit the cross sections reported by 
MiniBooNE [121 EH] , is implemented to linearly decrease 
to zero for the neutrino energy between 1 and 5 GeV, to 
avoid inconsistency with the results from NOMAD [5]. 
As the experimental data constrain the two-nucleon con¬ 
tribution for the carbon target only, for other nuclei its 
strength is assumed to exhibit linear dependence on the 
mass number A. 

In QE scattering, the default nuclear model in genie 
is the relativistic Eermi gas (REG) model of Bodek and 


Ritchie [50] . in which a high-momentum tail—inspired 
by the effects of nucleon-nucleon correlations—is added 
to the nucleon momentum distribution. 

Using the vT package [30], we replace the REG model 
by the spectral function (SF) approach [HI], the imple¬ 
mentation of which has been validated through compar¬ 
isons with electron-scattering data. In the SF formalism, 
the interaction between the beam particle and the nu¬ 
cleus is assumed to involve a single nucleon, with the 
remaining {A — 1) nucleons acting as spectators. Under 
this assumption—called the impulse approximation—the 
target initial state can be described by the nuclear spec¬ 
tral function P{p,E), giving the probability distribution 
that removal of a nucleon with momentum p from the 
nuclear ground state leaves the residual {A — l)-nucleon 
system with excitation energy E. 

The realistic SF of carbon, employed in the vT pack¬ 
age, has been obtained by the authors of Ref. m in the 
local-density approximation (LDA). The EDA scheme 
relies on the premise that surface and shell effects do 
not affect short-range correlations between nucleons in 
nuclei, and it is supported by the evidence that for 
IpI ^ 300 MeV the momentum distribution is largely 
independent of the mass number A, for nuclei with 
A > 4 |ST]. Therefore, using the LDA, the correlation 
contribution to the nuclear SF can be obtained from the¬ 
oretical calculations for uniform nuclear matter at differ¬ 
ent densities [3I1I52] and consistently combined with the 
shell structure of the nucleus, deduced from experimen¬ 
tal (e, e'p) data [53l iMj. We stress that the carbon SF 
of Ref. [3T] has proven successful in describing electron¬ 
scattering data in various kinematical setups, see e.g. 
Ref. [55] , and the momentum distribution obtained from 
it is consistent with the one extracted from (e, e'p) data 
at large missing energy and momentum [56] . 

The effect of using the lyT package on the total in¬ 
clusive GG cross section for muon neutrinos in genie is 
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FIG. 2. (color online). Same as in Fig. [^but for the muon antineutrino. 


presented in Fig. 1(a) On the one hand, the SF calcu¬ 


lation of the QE cross section improves the agreement 
with the experimental data at low energies reported by 
the SciBooNE Collaboration M- On the other hand, 
it does not alter the good agreement at the kinematics 
of the NOMAD experiment [8], where the inclusive cross 
section is dominated by the contributions of DIS pro¬ 
cesses and resonant pion production; see Fig. |l(b)[ Note 
that in the case of the NOMAD data, the error bars rep¬ 
resent A{a/E^), whereas for the SciBooNE points, they 
correspond to {Aa)/Ei, [Tlj . 

While the calculations do not seem to reproduce the 
SciBooNE point at the average neutrino energy of 2.47 
GeV, it is important to note that its energy range extends 
over the whole high-energy tail of the flux and does not 
end at 3.5 GeV. 


The recent results of the T2K experiment [iTl [19] sug¬ 
gest that in the ~1 GeV region uncertainties of the in¬ 
clusive cross section may be sizable, of the order of 20%, 
and that the contribution of the 2p2h processes may be 
smaller than that deduced from MiniBooNE. Should the 
2p2h estimate be significantly reduced in genie, the cal¬ 
culated cross sections would be in very good agreement 
with the T2K data. 

For the sake of completeness, in Fig.[^ we show the to¬ 
tal inclusive GG cross section for muon antineutrinos, for 
which experimental data are not available. It clearly ap¬ 
pears that the 2p2h contribution—which in genie m is 
assumed to be the same for neutrinos and antineutrinos— 
plays a much more important role in the latter case. 

As this article is focused on the effects of realistic detec¬ 
tion capabilities on energy reconstruction, we leave com¬ 
parisons to exclusive cross sections for future studies of 
nuclear-model uncertainties of the results presented here. 
Being able to reproduce the inclusive cross sections—of 
minimal uncertainties—the employed description of neu¬ 
trino interactions can be deemed suitable for the purpose 
of this analysis. 

In the context of the calorimetric energy reconstruc¬ 
tion, it is important to estimate how different species of 
hadrons contribute to the hnal-state energy. In Fig. 



FIG. 3. (color online). Fraction of the (a) and (b) 
energy converted into the kinetic energies of neutrons (solid 
lines) and protons (dotted lines) and the total energies of 
charged (short dashed lines) and neutral (long dashed lines) 
pions. 


we present the obtained fraction of the and energy 
carried out by neutrons, protons, and pions. 

While at energies below ~1 GeV the neutron results 
for neutrinos and antineutrinos clearly differ, this is not 
the case when DIS and resonant pion production become 
the dominant reaction mechanisms. Having in mind the 
numerical uncertainties—estimated not to exceed 1%— 
we conclude that at E^, > 1.5 GeV neutrons contribute 
less than 15% of the (anti)neutrino energy. 

As can be expected from QE interactions, the con¬ 
tribution of protons in neutrino interactions resembles 
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that of neutrons in antineutrino scattering and vice versa. 
For antineutrinos, the knocked-out protons carry out less 
than 10% of the initial energy in the whole considered re¬ 
gion. However, for neutrinos, the protons contribution is 
much more sizable and reduces to < 10% for energies 
exceeding ~3.35 GeV. 

Although the total energies of 7r°’s enter the calori¬ 
metric energy reconstruction, in the considered kinematic 
range, their contributions are smaller than that of neu¬ 
trons. On the other hand, in the region dominated by 
resonant and nonresonant pion production, charged pions 
give the largest contribution to the final-state hadronic 
energy. 

We observe that at energies above 0.5 GeV the analogi¬ 
cal results for electron neutrinos and antineutrinos do not 
differ significantly from those presented in Fig. 


IV. CONSIDERED DETECTOR EFFECTS 


In our analysis, two extreme sets of assumptions are 
considered regarding the detector performance: 

(i) Perfect reconstruction. 

With the exception of neutrons, all produced parti¬ 
cles are observed, and their measured energies are 
equal to the true ones. 

(ii) Realistic setup. 

The measured energies and angles are smeared with 
respect to the true ones by a finite detector resolu¬ 
tion. The detection efficiencies and thresholds are 
taken into account. Neutrons are assumed to es¬ 
cape detection. 

The detection thresholds applied in our calculations 
correspond to the measured kinetic energy of 20 MeV for 
mesons and 40 MeV for protons. For comparison, the 
NOMAD and MiniBooNE experiments were able to de¬ 
tect protons of kinetic energy above ~50 MeV [S], and the 
^40-MeV threshold is expected for future liquid-argon 
detectors m- 

The efficiencies are treated as energy independent, for 
the sake of simplicity, and set to 60% for 7r°’s, 80% for 
other mesons, and 50% for protons. Those values can be 
considered rather optimistic, compared to the efficiencies 
achieved in existing detectors [9l |58l |59] . We assume that 
produced charged leptons are always detected. 

Accounting for the effect of finite detector resolution 
on an observable ccmeas; we smear it according to the 
Gaussian distribution centered at the true value Xtme, 


y(^meas) — 


■\/2^fT(Xtrue) 


exp 


1 / ^meas :^true \ 

2 y (j(xtrue) / 


In the case of the muon, we apply this procedure to 
the momentum and production angle, using 


a(|k^|) =0.02|m and <j(0) = 0.7°, (12) 


as in the MINERvA experiment m- Eor the electron, 
the realistic resolutions 13 El] 

a(B^) = O.lOEe and a(0) = 2.8° (13) 


are used. 

For other particles, the only smeared quantity is the 
energy. To the electromagnetic showers produced by 
decays, we employ the energy resolution given by 


E^o 


= max 


OttO 1 

V%(o ’ Et,0 J 


(14) 


with = 0.107 and = 0.02, while for other hadrons. 


a(Eh) 

Eh 


= max 



(15) 


is used with au = 0.145 and bh — 0.067. The values of 
a and the energies appearing in Eqs. ( [I4| ) and (15) are 
expressed in units of GeV. Note that the hadron-energy 
resolutions applied in our analysis can be considered op¬ 
timistic, as they are to be compared to 


crjE^o) 

Et,0 



TT^ 

El, 


and 


(xjEh) 

Eh 



K 


achieved in the MINOS [4] and MINERvA experi¬ 
ments. Note that the assumed energy resolution is 
also optimistic when compared with the expectation for 
future liquid-argon detectors from Ref. [5^ . 

It is important to emphasize that in the context of the 
kinematic method of energy reconstruction our assump¬ 
tions are conservative. We employ realistic detector res¬ 
olutions for muons and make minimal use of information 
on observed hadrons. Disregarding additional protons 
in the final state, we reconstruct every event as if it in¬ 
volved single-nucleon knockout. All mesons in the final 
state are assumed to originate from a single resonance 
of the invariant mass Ma = 1.232 GeV, and their angu¬ 
lar distributions are not taken into account; compare to 
Eq. (|^ . On the other hand, detector effects entering the 
calorimetric method are treated in a rather optimistic 
way, compared to those in existing detectors. 

Note that knocked-out neutrons are assumed to be un¬ 
detected both in the realistic and perfect scenarios, owing 
to difficulties of their reconstruction in neutrino events. 
Traveling some distance from the primary interaction 
vertex before scattering, neutrons are currently problem¬ 
atic to associate with the neutrino event, and they typ¬ 
ically deposit only part of their energy in the detector. 
However, this assumption may need to be revisited in the 
future, when an ongoing experimental program [631164] 
brings important progress in the understanding of the 
detector response to neutrons. 


V. MIGRATION MATRICES 

In our calculations, detector effects are taken into ac¬ 
count in the migration matrices, AI^, the columns of 















7 



0.5 I.O 1.5 2.0 2.5 




0.5 1.0 1.5 2.0 2.5 



E,ec [GeV] 


Erec [GeV] 


FIG. 4. (color online). Probability distribution functions for a 
event of the energy iStrue to be reconstructed at an energy 
Erec- The results obtained for a QE event at i?true = 1-45 GeV 
(top) and a DIS event at iStrue = 3.45 GeV (bottom) show the 
effects of different assumptions on the detector performance; 
see the text for details. 


which are the probability distribution functions (PDFs) 
for a neutrino interaction X at the true energy in the fth 
bin to be reconstructed with an energy in the jth bin. 
The observed event distribution can then be calculated 
as 


V j 

where X runs over the four types of interactions consid¬ 
ered (DIS, res, 2p2h, and QE), i and j refer to the energy 
bins, and stands for the number of X events in the 
bin j computed without detector effects. Note that, bar¬ 
ring the boundary effects, the smearing produced by the 
migration matrices does not have any impact on the total 
number of events. 

To facilitate the reproduction of our results, we provide 
the complete set of migration matrices, calculated for the 
energies up to 8 GeV, using 0.1 GeV bins, and the cross 
sections employed in our analysis [55] . 

In an ideal detector, all migration matrices would be 
the unit matrices. In a real experiment, the reconstructed 
energy may depend on the reconstruction procedure, as 
well as on the mechanism of interaction. Imperfect de¬ 
tection capabilities—energy resolutions, efficiencies, and 


FIG. 5. (color online). Same as in Fig. |^but for a event. 


thresholds for particle detection—affect the probability 
for a neutrino event to be reconstructed in the correct 
energy bin. In particular, finite energy resolutions smear 
the measured energies, while imperfect efficiencies and fi¬ 
nite thresholds result in an energy partially carried away 
by undetected particles. As a consequence, PDFs have 
finite widths, they are asymmetric, with a broader tail 
toward the lower energies, and their mean values are low¬ 
ered with respect to the true neutrino energies. 


Those features are illustrated in Figs. and for the 
QE (DIS) mechanism of interaction and the true neutrino 
and antineutrino energy of 1.45 GeV (3.45 GeV) recon¬ 
structed using the calorimetric method. Gomparing the 
PDFs obtained for the perfect and realistic reconstruc¬ 
tions, defined in Sec. IV one can observe how sizable the 
effect of realistic detector capabilities is. For reference, 
we also show the Gaussian distribution with the standard 
deviation a{E^) = 0.15^/E^, with cr and in units of 
GeV, typically applied to account for detector effects in 
phenomenological studies devoted to liquid-argon detec¬ 
tors lilSTlIMiHj- 


The differences between the neutrino and antineutrino 
PDFs can be traced back to a twofold reason: (i) differ¬ 
ent contributions of neutrons to the final-state energy, as 
shown in Fig.[^ and (ii) the typical energy transfer being 
lower in antineutrino interactions, owing to the destruc¬ 
tive interference of the response functions. 

Even for the perfect-reconstruction scenario, the PDFs 
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FIG. 6. (color online). Effective energy resolution as a func¬ 
tion of the true energy for muon neutrinos. The results of our 
Monte Carlo simulations for QE (DIS) events are shown as 
lower (upper) bands. For each band, the upper (lower) edge 
corresponds to the realistic (perfect) detection capabilities, 
defined in Sec. HYl For comparison, a few smearing functions 
frequently used in phenomenological oscillation studies are 
also shown. 


are asymmetric, as a consequence of pion absorption in 
the nuclear medium and the energy carried out by neu¬ 
trons, assumed to escape detection. When the realis¬ 
tic detector effects are accounted for, the PDFs clearly 
broaden due to the employed energy resolutions, and 
their modes shift toward lower energies. The latter effect 
is particularly large for DIS events, in which the muon 
contribution to the final-state energy is typically smaller 
than in QE scattering, and the role of the efficiencies is 
larger. 

To make contact with existing phenomenological stud¬ 
ies of neutrino and antineutrino oscillations, in Figs. 
and we compare a few simple functions typically em¬ 
ployed as an effective energy resolution with our esti¬ 
mates based on the Monte Carlo simulations, for both the 
calorimetric (upper panels) and kinematic (lower panels) 
reconstruction methods. Our calculations of the PDF’s 
standard deviations are presented as bands spanning the 
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FIG. 7. (color online). Same as in Fig. [^but for 


values between the results for the perfect-reconstruction 
scenario (lower edge) and those for the realistic scenario 
(upper edge). The estimates for QE and DIS events 
are shown separately, as the lower (darker) and upper 
(lighter) bands. Because the energy resolutions applied 
in the realistic scenario exceed those in existing exper¬ 
iments roughly two times, they can be considered opti¬ 
mistic. On the other hand, an effective energy resolution 
higher than our result for the perfect-reconstruction sce¬ 
nario does not seem to be possible to achieve without the 
reconstruction of neutrons. 

Note that for liquid-argon experiments (such as DUNE 
or LBNO) or totally active scintillator detectors (e.g., 
NOz^A), the effective energy resolution is typically as¬ 
sumed to be cr{E^) = O{0.05-0.2)^/E^, in units of GeV; 
see, for instance, Refs. On the other hand, 

for experiments based on the Cherenkov technique (such 
as T2K), the resolution a{E^) = 0{0.08)y/E{^ is usually 
applied in phenomenological studies [571 ED]. 

Finally, in Figs. and we show how the modes 
of the reconstructed-energy distributions depend on the 
true value of energy, comparing the calorimetric and kine¬ 
matic reconstruction methods. The bands represent our 
Monte Carlo results for QE and DIS events, the lower 
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FIG. 8. (color online). Mode of the reconstructed-energy 
distributions as a function of the true energy calculated for 
muon neutrinos. The bands show our Monte Carlo results, 
with the lower (upper) edge obtained assuming the realistic 
(perfect-reconstruction) scenario. The darker (lighter) bands 
present the results for QE (DIS) events. For reference, the 
lines corresponding to the true value and its underestimation 
by 10 and 20% are also shown. 



E,„e[GeV] 
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FIG. 9. (color online). Same as in Fig. [^but for zJ^’s, 


(upper) edge of which corresponds the realistic (perfect- 
reconstruction) scenario. 

In the calorimetric method, the modes for QE scatter¬ 
ing are in a good agreement with the true energy, and 
the expected presence of neutrons in the final state of 
antineutrino events introduces only a small effect. While 
for DIS, the agreement is somewhat reduced—especially 
when detector effects are taken into account—the modes 
do not differ from the true energy by more than 20%. 

For the kinematic method, in QE interactions, both the 
neutrino and antineutrino modes are in excellent agree¬ 
ment with the true energy. However, this is not the case 
for DIS events. For antineutrinos, the discrepancy be¬ 
tween the DIS mode and the true energy exceeds 10% 
and is typically close to 20%. For neutrinos, the mode 
underestimates the true energy by ^1 GeV. This behav¬ 
ior can be traced back to the simplicity of our kinematic 
analysis, assuming that all events containing at least one 


pion are produced through excitation of a resonance of 
the invariant mass M\ = 1.232 GeV; see Eq. Should 
an experiment be able to reconstruct pions angles with 
high angular resolution, more refined treatment (|^ would 
be possible. 


The bands in Figs. [6]^ show the differences between 
the results for the realistic and perfect-reconstruction sce¬ 
narios. Their widths can be used as a measure of the sen¬ 
sitivity to detector effects, which turns out to be larger 
for the calorimetric migration matrices than for the kine¬ 
matic ones. However, they do not represent uncertain¬ 
ties entering any actual experiment. In practical situa¬ 
tions, detection capabilities and their uncertainties are 
estimated in test-beam exposures. In Sec. |VI1 we are go¬ 
ing to analyze the importance of an accurate estimate of 
these uncertainties in the context of an oscillation anal¬ 
ysis. 
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VI. OSCILLATION ANALYSIS 

The event number expected in a neutrino-energy bin, 
without considering detector effects, can be computed as 

Nl = / dE, a,,{E,) {{e},E,) 

J Ei 

(16) 

where the indices a and /3 label the initial and final neu¬ 
trino flavors, respectively; lS.Eu is the energy-bin size; 

denotes the vp cross section; and is the unoscil¬ 
lated Va flux. The Va —>■ vp oscillation probability, Pap, 
depends on a set of oscillation parameters { 0 } as well as 
on the neutrino energy E^,. In the oscillation 

channel, the oscillation probability can be approximated 
by 

~ 1 - sm 0^^ sm I 1 , 

where 


energies, either the accuracy is reduced [see Eq. ([^ con¬ 
taining unknown on event-by-event basis] or detailed 
angular information is required [see Eq. Qj. 

Phenomenological analyses for wide-band neutrino 
beams operating in the multi-GeV energy regime have 
shown that the sensitivity to neutrino oscillation param¬ 
eters of a liquid-argon detector is similar to that of a 
Cherenkov detector of the mass ^3-6 times larger, be¬ 
cause of the lower efficiency of the latter for non-QE 
events [571 UnHHS]- It has been argued that the rea¬ 
son lies in the imaging capabilities of the liquid-argon 
technology, able to identify protons in addition to lepton 
tracks m, as opposed to Cherenkov detectors in which 
the information on CC events comes predominantly from 
the charged leptons. Nevertheless, in such comparisons, 
the neutrino energy resolution at liquid-argon detectors 
was always assumed to be extremely good, in the range 
of a{Ey) = QAby/Edy. Should this be affected by detec¬ 
tor effects, or by the energy carried away by unobserved 
particles, such conclusions may have to be reexamined. 


sin^ = 4cos^ 6 »i 3 sin^ 6»23(1 - cos^ di 3 sin^ 6 * 23 ), 

and is the weighted average of Am|]^ and 

Am '^2 |H] that can be expressed as 

Aml^ = Amli + 0{Aml^). 

Should the experimental information on the Ey distri¬ 
bution of events be available, the values of the oscillation 
parameters {CH could be, in principle, directly inferred 
from the datrR However, modern neutrino beams are 
produced as tertiary products, originating predominantly 
from the decay of pions and kaons produced in interac¬ 
tions of protons impinging on a target. Therefore, neutri¬ 
nos are not monoenergetic, and—for a given event—the 
incident neutrino energy has to be reconstructed from 
the kinematics of the particles in the final state. Sys¬ 
tematic uncertainties of this procedure inevitably depend 
on the capabilities of the employed detector, as well as 
on the neutrino interaction channel, owing to the differ¬ 
ent detection efficiencies for the particles involved in the 
event. Because for muons the efficiency uncertainty is 
minimal, its effect on the reconstructed energy spectrum 
is neglected in our analysis. 

The kinematic method of energy reconstruction—used, 
for instance, in Cherenkov detectors—is known to be ac¬ 
curate to ~I00-I50 MeV for QE events (T^. However, 
as discussed in Sec. |H| this is not the case for the events 
of QE topology containing undetected hadrons. For ex¬ 
ample, a single undetected pion typically spoils the re¬ 
constructed energy by ~300-350 MeV [35]. For pion- 
production events, playing an important role at higher 


^ Note, however, that if more than one parameter is being deter¬ 
mined from the data, severe degeneracies among the different 
parameters may take place. I72H78I 


In this work, we compare the kinematic and calori¬ 
metric methods of neutrino energy reconstruction in the 
context of a disappearance experiment. Four types of 
CC neutrino interactions are considered—resonant and 
nonresonant pion production, two-nucleon knockout, and 
quasielastic scattering—and modeled according to genie 
with the eT package, as described in Sec. |TTT| The coher¬ 
ent channel has not been taken into account because of its 
negligible contribution to the event rate [84]. For each 
interaction type and both reconstruction methods, we 
have calculated corresponding migration matrices 
with and without detector effects; see Secs. IV and[^ 


In the —>■ oscillation channel, the main back¬ 

ground comes from neutral current events misidentified 
as CC ones. As it is generally expected to be very low, 
for simplicity, we neglect it in our analysis. On the other 
hand, due to the large event statistics, systematic un¬ 
certainties are expected to be relevant. In their treat¬ 
ment, we follow that of Refs. [551135], with the imple¬ 
mentation as detailed in the appendices of Refs. [551157] . 
A 20% bin-to-bin uncorrelated systematic uncertainty is 
assumed, as well as a 20 % overall normalization uncer¬ 
tainty, which is bin-to-bin correlated. The pull method is 
used, adding a Caussian prior for each systematic error, 
and the final profile is obtained after the minimiza¬ 
tion over the nuisance parameters. These (rather con¬ 
servative) systematic uncertainties are introduced in our 
analysis in order to accommodate possible differences in 
the shape of the expected event distributions at the de¬ 
tectors, in a similar fashion as was done in Refs. [851l86j . 
A near detector is also considered in the analysis, which 
helps to constrain the nuisance parameters during the fit. 


The oscillation analysis is performed using a modified 
version m of GLOBES [551 inn] ■ The assumed true values 
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of the oscillation parameters are 

= 7.50 X 10"® eV^; = 2.46 x lO'^ eV^; 

012 = 33.48°; 023 = 42.30°; 0x3 = 8.50°; (5 cp = 0. 

(17) 

In the analysis, we focus on the determination of at¬ 
mospheric parameters 023 and Arri^i- For simplicity, all 
oscillation parameters have been fixed during the anal¬ 
ysis; i.e., no marginalization has been performed when 
obtaining the allowed confidence regions. However, our 
conclusions are not expected to change significantly if 
marginalization was performed within the currently al¬ 
lowed experimental regions. Matter effects are included 
in our simulations, and the matter density profile has 
been chosen according to the preliminary reference Earth 
model [90] . 

A. Considered experimental setups 

Oscillation experiments using neutrino beams pro¬ 
duced from meson decays in flight can be divided into 
two main categories according to their far detector loca¬ 
tions: on-axis (such as K2K and MINOS [3] and also 
the DUNE [6| and LAGUNA-LBNO |92] proposals) and 
off-axis (for instance, ongoing T2K [T] and NOpA IS])- 

Because of the pion-decay kinematics, the flux at an 
off-axis site is well localized around a given neutrino en¬ 
ergy, with the beam spread and high-energy tail heavily 
reduced [93]. Such a design allows for a signihcant re¬ 
duction of the backgrounds coming from neutral-current 
events, and guarantees that the range of values of L/E^, 
is well localized around the first oscillation maximum. 
The price to pay is a lower beam intensity compared to 
an on-axis configuration. In addition, as the mean neu¬ 
trino energy is lower than in an on-axis experiment, the 
number of events in the far detector is also lower, since 
the cross section is an increasing function of the energy. 

On the other hand, in an on-axis neutrino oscillation 
experiment, the spread of the beam is larger. In principle, 
this allows one to determine the shape of the oscillation 
probability by performing measurements at different val¬ 
ues oi L/E^. In addition, thanks to the higher beam 
intensity and higher neutrino energies, large event statis¬ 
tics is easier to collect. As a consequence of being more 
energetic, though, the event sample generally contains a 
significant fraction of pion-production events. In addi¬ 
tion, the high-energy tail typically produces a significant 
neutral-current background contaminating in particular 
the low-energy bins. 

In this article, we perform an analysis of two differ¬ 
ent neutrino-oscillation experiments: a HE setup with a 
broadband on-axis beam, and a LE setup with a narrow- 
band off-axis beam. 

In the LE setup, we use the neutrino beam of Ref. [94], 
being a preliminary estimate of that in the T2K experi¬ 
ment with the 2.5° off-axis configuration [53]. As shown 



E„[GeV] 

FIG. 10. (color online). The muon neutrino fluxes (in arbi¬ 
trary units) as a function of the neutrino energy for the two 
configurations considered in this work. The solid (dashed) 
line corresponds to the low-(high-)energy setup. 


TABLE I. Expected number of unoscillated CC events for 
the low- and high-energy setups defined in the text. In addi¬ 
tion to the total event numbers, their breakup into different 
reaction mechanisms, labeled as in Fig. [^ is also given. 



QE 

2p2h 

res 

DIS 

total 

LE (0.3-2 GeV) 

49% 

28% 

21% 

2% 

4891 

HE (0.3-4 GeV) 

26% 

11% 

37% 

26% 

4456 


in Fig. [^ it is peaked at around 600 MeV. For the 
distance to the far detector for the LE setup, we take 
295 km, the length corresponding to the first oscillation 
maximum. The normalization (i.e., number of protons 
on the target, detector size, and data-taking time) is ar¬ 
bitrarily set to obtain a total unoscillated CC inclusive 
event sample of 0(5000) events in the energy range be¬ 
tween 0.3 and 2 CeV, relevant to the analysis. 

For the HE option, the neutrino flux labeled as 
“650 km” in Ref. [55] is used, being one of the configura¬ 
tions considered within the LACUNA design study m- 
This flux has a broad peak between ~1 and ~2 CeV, 
with a non-negligible tail extending well above 3 CeV; see 
Fig. [^ We considered several choices for the distance 
to the far detector in the range L = 0(500-1000) km, 
for which the hrst oscillation maximum would lie within 
the energy range of the flux peak, and found that L ^ 
1000 km gives optimal results in the disappearance 
channel. In the following, we discuss our results for this 
baseline choice only. 

Since in this work we are only interested in exploring 
detector effects on different energy-reconstruction meth¬ 
ods, the normalization for the HE setup is arbitrarily set 
to obtain a similar total number of unoscillated events 
as for the LE setup. However, due to the much higher 
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neutrino energies, the composition of the inclusive event 
sample is very different. To illustrate it, the total number 
of unoscillated CC inclusive events at the far detector are 
given in Table for both setups, together with the per¬ 
centages for QE, 2p2h, res and DIS events in each sample. 

In our analysis, all neutrino events with the energies 
between 0.2 and 8 GeV are calculated for both LE and 
HE setups. However, only those which are reconstructed 
between 0.3 and 2 GeV (4 GeV) for the LE (HE) setup are 
considered during the fit. The is built by binning the 
events in reconstructed neutrino energy, using 100 MeV 
bins. 


B. Results for the calorimetric method 


As explained in Sec. |TTj in the calorimetric method, 
the neutrino energy is reconstructed by adding the en¬ 
ergies of all observed particles in the final state, and no 
information on the direction of the outgoing particles is 
used. 

In an analysis of an oscillation experiment, one has 
to rely on a Monte Garlo simulation to predict the “ex¬ 
pected” (or “fitted”) event rates. One of the inputs to 
such a simulation is detailed information on the detector 
capabilities. However, it is subject to uncertainties, the 
role of which we are going to analyze. 

Our results for the calorimetric method are shown in 
Figs. E] and [T^ for the LE and HE setups, respectively. 
In their upper panels, we present the simulated event 
distributions in the far detector as a function of the re¬ 
constructed energy. The lower panels illustrate the con¬ 
fidence regions for fits to the atmospheric oscillation pa¬ 
rameters. In the fits, the true event rates are computed 


for the oscillation parameters (17) taking into account 
detector effects according to the realistic setup described 
in Sec. ITVl 

Because it is computationally expensive to generate a 
given set of migration matrices from simulated neutrino 
events, in our analysis, the fitted rates are generated us¬ 
ing a linear combination of the matrices obtained for the 
realistic and perfect-reconstruction scenarios. Allowing 
for a continuous transformation of the event-rate distri¬ 
bution from one scenario to the other, this method is 
useful to quantify when the incorrect estimation of the 
detector performance starts to have a significant impact 
on the fit. In practice, the fitted rates are obtained as: 


= 


5] ^ {(1 - a)M^' } TVf, 

^ .7 


( 18 ) 

where X stands for a given type of interaction (QE, 
2p2h, res, or DIS) and i and j are the energy-bin indices. 
]\px, real j^jX, perf the migration matrices ob¬ 

tained for the realistic and perfect-reconstruction scenar¬ 
ios, respectively. Finally, is the event rate in the en¬ 
ergy bin j for the interaction type X, computed without 
accounting for the detector effects, as in Eq. (16). 
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FIG. 11. (color online). Resnlts for the low-energy setnp and 
the calorimetric reconstruction method. Upper panel: Simu¬ 
lated charged-current event distributions in the far detector 
as a function of the reconstructed energy, for the oscillation 
parameters (17l. The shaded histogram corresponds to a per¬ 
fect estimate of detector effects (see Sec. IV I. The dashed and 
solid lines show the event rates obtained with the detector 
performance overestimated by 10% and 30%. Lower panel: 
Confidence regions in the (023, Amii) plane, at Icr C.L. (2 de¬ 
grees of freedom). The shaded area corresponds to the perfect 
estimate of the detector effects. The lines show the contours 
obtained for the detector performance overestimated by 10%, 
20%, and 30% in the fit. The star indicates the true values of 
the oscillation parameters, the same for all confidence regions 
shown. 


In Eq. (18), a is a purely phenomenological parameter, 
used to obtain the “effective” migration matrices as a 
linear deformation of the two extreme scenarios under 
consideration. For instance, a = 0 means that the fitted 
rates are obtained in the same way as the true rates. On 
the other hand, a = 0.3 means that the fitted rates are 
obtained underestimating the role of detector effects by 
30%, a substantial amount. Two examples of the event 
distributions obtained for different values of a are shown 
in the upper panels in Figs. [TT] and [T^ Even though the 
distributions look very different, all histograms have been 
obtained assuming the same true values for the oscillation 
parameters ( [T7| . 

For a given value of a, the fitted rates are simulated for 
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FIG. 12. (color online). Same as in Fig. 11 but for the high- FIG. 13. (color online). Same as in Fig. 11 but using the 


energy setup and the calorimetric reconstruction method. 


kinematic method to reconstruct the neutrino energy. 


every possible combination of (Am|]^, 023). The point in 
the plane which gives a best fit to the data 

is then identified, and the confidence regions are drawn 
by requiring that 


Ax (Ato3i, 6»23) = X (Am3i,6l23) - Xbest-flt < 2-30. 

(19) 

Since the fitted rates are generated using a different set of 
the migration matrices, the event distributions generally 
have a different shape than the true rates, and the best 
fit to the data does not necessarily coincide with the true 
values of oscillation parameters realized in nature. When 
the detector capabilities are incorrectly estimated, the 
allowed confidence regions start to drift away from the 
shaded areas, as shown in the lower panels in Figs. 
and 12 The minimum value of the x^ obtained for each 
a is indicated in the legend, together with the effective 
number of degrees of freedom (d.o.f.) in the fit, Ad.o.fj^ 
For a = 0 the value of the x^ at the best-ht point is 


consistent with zero and, therefore, is not indicated in 
the legend. 


C. Results for the kinematic method 


In the kinematic method, the neutrino energy is recon¬ 
structed using only the direction and momentum of the 
produced charged lepton in the final state, and no kine¬ 
matic information from the outgoing hadrons is needed. 

As in Sec. VI B| we calculate the true event distribu¬ 
tion using the migration matrices accounting for realis¬ 
tic detection capabilities, while the fitted rates are ob¬ 
tained following Eq. (18), and the confidence regions sat¬ 
isfy Eq. ([T^. 


Our results for the kinematic reconstruction method 


are shown in Figs. and 14 for the LE and HE setups. 


^ The effective number of degrees of freedom, being the number 
of bins used in the fit minus the number of parameters deter¬ 
mined from the data, is typically used in tests of the validity of 
the model used to fit the data. In a real experiment, the mini¬ 
mum would receive contributions from statistical fluctuations 


of the measured quantities in each data bin, and XminZ-^d.o.f. 
should be close to 1, if the model used to fit the data is correct. 
In our case, no statistical fluctuations are considered in the fit. 
Nevertheless, the use of a wrong model in data fitting can give 
a sizeable contribution to the value of Xmin/^d.o.f. • This can be 
estimated in our calculations as Xbest-fit/'^d.o.f.- 
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FIG. 14. (color online). Same as in Fig. 11 but for the high- 


energy setup and the kinematic energy reconstruction. 


respectively. In the case of the confidence regions rep¬ 
resented by the shaded areas, the true and fitted rates 
are generated using the same set of migration matrices 
{a = 0), and, therefore, the best fit coincides with the 
true values. When the value of a increases, the confi¬ 
dence regions start to drift away from the shaded areas. 
This is a consequence of the fact that the shapes of the 
event distributions are generally different from that of 
the true rates, even for the same oscillation probability, 
and the best fit for the oscillation parameters no longer 
coincides with the true values. 

However, it should be noted that for the kinematic 
reconstruction, the effect of the mismatch between the 
actual and assumed detector capabilities is much milder 
than for the calorimetric method. The confidence regions 
show a significant overlap, even for a = 0.3, correspond¬ 
ing to the detector performance substantially overesti¬ 
mated. The main reason for this is that the muon re¬ 
construction is very precise in modern detectors, while 
for electrons and hadrons this is more difficult to achieve 
(see also Figs. and §. It should be noted, however, 
that in this study we assume that the muon track is fully 
contained in the detector, regardless of its energy. Should 
this assumption be relaxed, the muon momentum deter¬ 
mination may be affected and, as a consequence, also the 
neutrino-energy reconstruction. 


VII. SUMMARY 

It has been realized for a while now that an accurate 
understanding of neutrino-nucleus scattering is an essen¬ 
tial ingredient for neutrino-oscillation experiments em¬ 
ploying accelerator beams. More recently, quantitative 
studies have tackled the relation between uncertainties 
in cross section modeling and the resulting physics sensi¬ 
tivities for oscillation measurements. For example, it was 
shown that for quasielastic events the energy reconstruc¬ 
tion based exclusively on the kinematic observables of the 
outgoing lepton is susceptible to a large bias resulting 
from the underlying nuclear-interaction model [851186] . 
In this article, we fully rely on the nuclear model used 
for event generation, and have a critical look only at the 
impact of realistic detector effects. 

From energy conservation, it is clear that a perfect 
calorimeter—i.e., a detector able to measure the total 
energy of all reaction products—would be free from any 
bias in energy reconstruction [98j . In this article, we 
aim at quantifying how close to the perfect scenario the 
calorimetric energy reconstruction is, when finite energy 
resolutions, detection efficiencies, and thresholds are ac¬ 
counted for. In addition, we compare the role of realistic 
detection capabilities on the calorimetric and kinematic 
analysis. In the latter case, we use—although in a sim¬ 
plistic way—the information on that the charged lepton’s 
kinematics also for non-QE events. 

We find that the kinematic reconstruction is very ro¬ 
bust with respect to detector effects, largely because 
muons are the most precisely reconstructed particles in 
modern neutrino detectors. On the other hand, the ac¬ 
tual performance of the calorimetric analysis clearly de¬ 
pends on the assumed detector performance. The em¬ 
ployed detection capabilities are not meant to represent 
any existing detector, but are indicative of the general 
level of performance which can be expected. While to 
translate these result into a specific experiment a de¬ 
tailed study of detector response—beyond the scope of 
this work—would be required, many of the analysis tech¬ 
niques developed here should turn out to be very useful. 
Interestingly, we find that the kinematic reconstruction 
performs well even for pion-production events, but the 
independence of this observation from the underlying nu¬ 
clear model is to be examined. 

We limit our study to the disappearance channel 
since this can be effectively treated as a two-flavor os¬ 
cillation. As a consequence, the effects of energy reso¬ 
lution and misreconstruction are directly related to the 
precision of the Am^ determination and a shift of the 
value of Ato^, respectively. We use a phenomenologi¬ 
cal parametrization to interpolate between a perfect and 
realistic detector. This allows us to conclude that, over¬ 
all, the detector response—in terms of efficiencies, resolu¬ 
tions, and thresholds for individual particles—has to be 
understood at a 10% level or better, to avoid a significant 
bias in the measurement of Am^. For the kinematic anal¬ 
ysis, these requirements are much less stringent, but un- 
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certainties of the nuclear model—not considered here— 
present a challenge, as shown previously in the literature. 

In summary, while the calorimetric reconstruction may 
be less sensitive to the underlying nuclear model, it is 
strongly affected by detector effects, typically leading to 
energy underestimation. On the other hand, the kine¬ 
matic method of neutrino energy reconstruction is much 
less challenging for the detector design, but it strongly 
relies on an accurate understanding of neutrino-nucleus 
interactions. One needs to keep in mind that the quan¬ 
titative details of our conclusions may be specific to 
this work, owing to the underlying detector assumptions. 
However, their qualitative aspects can be expected to 
hold for a variety of experiments. 

As a final remark, we observe that the results presented 
in this article are subject to nuclear-model uncertainties, 
likely to be largest for 2p2h processes. Their estimate is 
left to be quantified in our future studies. 
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